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Abstract. We apply a mass reconstruction technique to 
large-scale structure gravitational distortion maps, simu- 
lated for different cosmological scenarii on scales from 2.5 
arcmin to 10 degrees. The projected mass is reconstructed 
using a non-parametric least square method involving the 
reduced shear on which noise due to intrinsic galaxy el- 
lipticities has been added. The distortion of the galax- 
ies is calculated using the full lens equation, without any 
hypothesis like the weak lensing approximation, or other 
linearization. 

It is shown that the noise in the reconstructed maps 
is perfectly uncorrelated Poissonian, with no propagation 
from short to large scales. The measured power spectrum 
and first four moments of the convergence can be corrected 
accurately for this source of noise. The cosmic variance 
of these quantities is then analyzed with respect to the 
density of the background galaxies using 60 realizations 
of each model. We show that a moderately deep weak 
lensing survey (5x5 degrees with a typical background 
population of 30 gal/arcmin 2 at a redshift z s ~ 1) is able 
to probe the amplitude of the power spectrum with a few 
percent accuracy for models with a% Q°' s = 0.6. 

Remarkably, we have found that, using the third mo- 
ment of the local convergence only, such a survey would 
lead to a 6 a separation between open (f2 = 0.3) and flat 
(Q = 1) models. This separation does not require a very 
deep survey, and it is shown to be robust against different 
hypothesis for the normalization or the shape of the power 
spectrum. 

Finally, the observational strategy for an optimal mea- 
surement of the power spectrum and the moments of the 
convergence is discussed. 
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1. Introduction 

Mass reconstruction from gravitational distortion inver- 
sion is a promising technique to probe the mass distribu- 
tion and the clustering on very large scales, regardless of 
the nature and the dynamical state of the dark and lumi- 
nous matter. Pioneering theoretical work done by Gunn 
(1967), Jaroszyhski et al. (1990), Blandford et al. (1991), 
Miralda-Escude (1991) and Kaiser (1992) has shown that 
the expected distortion amplitude of weak lensing effects 
produced by large scale mass fluctuations (> IMpc) is 
roughly at the percent level. This low level of distortion 
is observable due to the large number density of galaxies 
observed in deep surveys (Kaiser 1992). Needless is to say 
that the observation of such distortion fields provides a 
unique way to build a picture of the large scale mass dis- 
tribution, independent on any biasing and/or dynamical 
prescriptions. 

The scientific impact on the determination of cosmo- 
logical parameters from weak lensing surveys has been un- 
derlined by recent theoretical work. Most of the papers 
quoted above show that the distortion two-point corre- 
lation function can be used to constrain the mass power 
spectrum. Villumsen (1996) remarked that the amplitude 
of the local distortion is proportional to the amplitude 
of the 3D density fluctuations and roughly to the density 
parameter £lo- Bernardeau et al. (1997, hereafter BvWM) 
extended these calculations in the S7o — A plane and show 
that the amplitude is also slightly dependent on the cos- 
mological parameter A. In the same paper the authors also 
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show that the shape of the convergence probability dis- 
tribution function can be used to disentangle the Q and 
as dependence for models of large-scale structure forma- 
tion with Gaussian initial conditions. In particular they 
demonstrate that the skewness 1 , third moment expressed 
in terms of the square of the second, is roughly inversely 
proportional to the total mass density of the Universe Qq, 
but independent on the amplitude of the fluctuations, as 
well as the shape of the power spectrum. This result, ob- 
tained by means of perturbation theory, is expected to be 
exact at large enough scale (see Gaztahaga & Bernardeau, 
1998). The skewness should be enhanced in the nonlinear 
regime at small scales (Gaztahaga & Bernardeau 1998, 
see also Colombi et al. 1997, and Jain et al. 1998) which 
amplifies the differences between open and flat cosmolo- 
gies. So, even if definitive quantitative predictions for such 
a quantity cannot be given from our present knowledge, 
it is clear that the skewness can accurately discriminate 
between different cosmological models. 

Recently many theoretical and observational aspects 
have been investigated in detail in order to converge to- 
wards an unbiased measurement of such small distortions 
(Bonnet & Mellier 1995, Kaiser et al. 1995, and Van Waer- 
beke et al. 1997). The technical limitations that have been 
recognized so far come from systematics due to the correc- 
tion of spurious distortions generated by the optic defects, 
fuzzy-shaped point spread function (PSF), and pixel con- 
volution (sampling). Available image analysis techniques 
and image quality ensure that such systematics can now 
be reduced to a one percent level, which is necessary for 
these weak lensing applications, but beyond the need for 
cluster mass reconstruction. The ultimate limitation for 
the use of weak lensing surveys as a cosmological probe 
depends therefore on the accuracy with which the correc- 
tion of the spurious distortions can be corrected. On the 
other hand, if the systematics can be reduced to the sub- 
percent level, the weak lensing analyses are limited by the 
intrinsic ellipticities of galaxies which acts as a shot noise 
for the gravitational distortion effect. 

In this paper we investigate the effects of the shot noise 
and finite size survey (cosmic variance) on the determina- 
tion of the cosmological parameters and the power spec- 
trum. In particular, two goals are looked for, 

— the precision on the determination of some cosmologi- 
cal parameters that can be expected from such maps; 

— the maximum level of observational systematics which 
is acceptable in order to achieve these theoretical pre- 
cisions. 

This study is made in particular in the perspective of 
ongoing and future wide field deep imaging surveys de- 
voted to weak lensing analysis (for instance the MEGA- 
CAM project, Boulade et al. 1998) at CFHT 2 or the SDSS 



project, Stebbins 1996). In such a perspective, there are 
still open issues concerning the optimal observing strategy 
that should be adopted: 

— Wide and shallow survey versus deep and narrow. 

— The determination of the optimal shape and size of 
the survey depends on the noise properties, and on the 
correlation properties of the signal. The investigations 
that have been done so far (Kaiser 1992, 1998 and 
Seljak 1997) assume that the projected mass follows a 
Gaussian statistics. If such an assumption had to be 
dropped it may significantly change the conclusions 
that have been reached. 

In addition, since the reconstruction of the projected 
mass from the shape of the galaxies is neither local nor lin- 
ear in terms of the distortion field (Kaiser 1995) and of the 
intrinsic ellipticities of the galaxies, it is essential to under- 
stand how the noise propagates in the reconstructed mass 
maps. In order to investigate how these different effects 
may couple together, we built a series of projected mass 
maps that contain a realistic amount of non-Gaussianity. 
The associated distortion field is then derived on the ba- 
sis of the full non-linear lens equation. A noise is added 
on the distortion maps, and the convergence is finally re- 
constructed. This paper presents the statistical analysis of 
those reconstructed maps for different cosmological mod- 
els and different observational contexts. 

The details of mass map generation and useful def- 
initions arc presented in Sect. 2. Sect. 3 presents the 
power spectrum analysis, the noise properties in the re- 
constructed mass maps, and the cosmic variance on the 
estimated power spectrum. Sect. 4 repeats similar analysis 
on moments in real space, where a comparison between top 
hat and compensated filter is done. It contains also some 
highlights about other possible statistical quantities that 
could be used to measure fi, and a comparison with results 
obtained for different power spectra. We finally summarize 
our results and discuss the best observational strategies. 

2. Generation of realistic K-maps 

2.1. Lensing effects, displacement and amplification ma- 
trix 

Any mass concentration deflects light beams by an an- 
gle proportional to the gradient of the local gravitational 
potential. This effect induces an apparent displacement 
of the sources, so that a source that was at the angular 

S I 

position £ will be observed at position £ , with 



2 V K (Xs-x) 
c 2 V K ( X )V K (xs) 



JdxWfat 1 ), (!) 



1 results are explicitly given at beginning of Sect. 4 

2 See also the WEB 
http://terapix.iap.fr/terapix_megacam.html. 
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where \ i s the radial coordinate (\s 1S the one of the 
source), T>k(x) t ne co- moving angular diameter distance 
and <p(x, £) the 3D gravitational potential. The differential 
displacement of the images induces an image distortion, 
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which depends on the second derivatives of the gravita- 
tional potential, i.e. on the mass density and the tidal field. 
This gravitational lensing effect of a thin lens is therefore 
characterized by an isotropic stretching, described by the 
convergence k, and an anisotropic distortion given by the 
complex shear 7 = 71 + 272. The so-called amplification 
matrix A describes the change of local coordinates be- 
tween the source and image planes and can be written 
as 



1 



K - 71 

-72 



1 



-72 

K + 7l 



(2) 



Its elements are related to the first derivative of the dis- 
placement field, 



71 = \ (V 2 £ 2 S - Vi£f ) 

72 = -Vi£f, 



(3) 

(4) 
(5) 



and are therefore related to the projected gravitational 
potential of the lens ip via, 



K = Aip/2; 71 = (ip n - V22)/2; 72 = ^12- 



(6) 



and where tp is given by, 
Jo 



where ip is the angle such that £ = T>k(x) 1 P- Here ip de- 
pends on x.s since the potential is integrated from a given 
source plane Xs to the observer. This can be trivially gen- 
eralized to a source redshift distribution. This equation 
is valid only if the lens-lens coupling is dropped and the 
Born approximation is used. Th validity of this assump- 
tion has been discussed by BvWM and Schneider et al. 
1997 (hereafter SvWJK). 

2.2. The galaxy shape matrices 

A source galaxy may be described by a complex ellipticity 
defined as 



l + r 



(8) 



where 1? is the galaxy orientation, and r is the square root 
of the ratio of the eigenvalues of the shape matrix Q*f) of 
the (centered) surface brightness profile of the galaxy, 



Q (s) = 



d0 9,6 3 S{0)/ J dOS(0). 



(9) 



In presence of lensing, the shape matrix of the image of the 
galaxy is given by Q — A~ 1 Q^A~ 1 , and provided that 
the amplification matrix does not vary over the galaxy's 



area, the observed ellipticity is still described by Eq. (9), 
and it writes (Schneider & Seitz 1995), 



e = 



+ 9 



1 + g*e( s ) ' 



(10) 



where g = 7/(1 — k) is the complex reduced shear. The 
observable is the distortion 5 = 2g/(l + \g\ 2 ), but for sub- 
critical lenses like those we discuss in this work, g is also 
directly observable, 



(ii) 



If the orientation of the source galaxies is random ((e^) = 
0), the observed mean ellipticity of galaxies is an unbiased 
estimate of the reduced shear (Schramm & Kayser 1995, 
Schneider & Seitz 1995), 



(e}=9. 



(12) 



The weak lensing approximation ((e) ~ 7) is generally 
used in the case of lensing by large scale structures. How- 
ever, we will not use that since we want to analyze the 
noise propagation in mass map reconstructed in the full 
non-linear regime (Eq. (10)). 

2.3. Construction of the synthetic projected mass maps, 
and their reconstruction 

The adopted procedure to build and analyze projected 
density maps is: 

1) Generation of a density map, which directly pro- 
vides the convergence k given the source redshifts: 

In order to make a precise analysis of the cosmic vari- 
ance and noise properties in the reconstructed mass maps 
it is necessary to have a large set of simulations of large 
scale structure. Instead of CPU intensive iV-body codes 
we use numerical fast second order Lagrangian dynamics 
(Moutarde et al. 1992) which has been shown to accu- 
rately reproduce the statistical properties of LSS (Munshi 
et al. 1994, Bernardeau et al. 1994, Bouchet et al. 1995). 
This allow us to build a large amount of projected mass 
maps for different cosmological models. The appendix A 
gives details about the generation of these maps and com- 
pares the non-Gaussian features with expectations from 
real dynamics. 

2) Calculation of the associated reduced shear map g 
and addition of a given level of noise to the reduced shear 
map: 

From the previous projected mass maps, a gravita- 
tional distortion map is computed, and the noise due to 
intrinsic ellipticities is added, as described in Appendix 
B. A "high" and "low" noise levels are simultaneously 
considered which correspond respectively to a mean num- 
ber density of galaxies of n = 30 gal/arcmin 2 and n = 
50 gal/arcmin 2 . For the I-band, these number densities are 



4 



Efficiency of weak lensing surveys 



reachable with respectively 1.5 or 4.5 hours exposure at 
CFHT which is expected to provide galaxies of redshift of 
about unity. At this stage we have at our disposal a large 
number of maps that are supposed to mimic accurately 
what can be observed with large CCD cameras. 
3) Reconstruction of the original k map from e: 
Many mass reconstruction methods already exists 
(Seitz ct al. 1998, Squires & Kaiser, 1996, and refer- 
ences therein). In this work, we use a non-parametric least 
square method (see Bartelmann et al. 1996 for details) 
which overfits data if no rcgularization process is used. 
The hope by doing so is to preserve all the noise proper- 
ties intact, so that a detailed noise analysis can be done. 
The appendix B gives more technical details about the 
reconstruction algorithm. 

Most of our generated images have an angular size of 
5x5 degrees, 120 x 120 pixels, each pixel having a 2.5' 
angular side (hereafter, the superpixel size) 3 . These are 
typical scales that a MEGACAM survey at CFHT could 
probe. Images of size 10 x 10 degrees (240 x 240 pixels) 
have also been generated, in order to estimate how the 
cosmic variance depends on the survey size in the non- 
linear regime, and to compare the merit of deep-small sur- 
vey area versus shallow-large survey. For each cosmologi- 
cal model and observational context series of 60 compact 
maps are produced. Table 1 summarizes the cases that 
investigated. 



Table 1. List of simulations which have been carried out. The 
power spectrum BG corresponds to the formula (38). It is the 
same for flo — 1 and Qo = 0.3. The CDM model corresponds 
to a standard CDM with T = 0.5. U is the ratio between the lo- 
cal convergence and the projected normalized density contrast, 
k = ui &2D (see appendix A). 
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•^sources 




size (deg 2 ) 


BG 


1.0 


0.6 


1.0 


0.115 


5x5 


BG 


0.3 


1.53 


1.0 


0.045 


5x5 


BG 


1.0 


0.6 


1.5 


0.195 


5x5 


BG 


0.3 


1.40 


1.5 


0.0837 


5x5 


BG 


1.0 


0.6 


1.0 


0.115 


10 x 10 


BG 


0.3 


1.53 


1.0 


0.045 


10 x 10 


BG 


1.0 


1.0 


1.0 


0.115 


5x5 


sCDM 


1.0 


0.6 


1.0 


0.115 


5x5 



Fig. 1 shows two examples of initial k maps, and the re- 
constructed mass of the noisy distortion maps. This panel 
illustrates what MEGACAM should be able to get dur- 
ing only 5 nights (!): 25 exposures in the I band, 1.5 hour 
each. A single MEGACAM field corresponds to a size of 

3 The reason for that pixelisation is that second order La- 
grangian dynamics is not able to probe high density peaks (i.e. 
small fields. 



24 x 24 pixels on Fig. 1, which is clearly the required mini- 
mum area to detect large scale structure features like super 
clusters, filamentary structures or voids. 

3. Power spectrum analysis of the reconstructed 
maps 

3.1. Noise statistical properties of the reconstructed maps 

It is clear that the mass reconstruction process does not 
produce any boundary effects (which is settled by defini- 
tion). The only boundary effect, slightly detectable on the 
figure, is a larger level of noise at the edge of the field, 
due to the change in the finite difference scheme at that 
position (see Appendix B). The noise due to the intrinsic 
ellipticities of the galaxies is clearly visible at small scales. 

Since the least x 2 method used to reconstruct the con- 
vergence is a local process, it is unlikely that noise prop- 
agates on scales larger than the pixel size 4 . 

Fig. 3 shows the power spectrum analysis of 60 recon- 
structed mass maps in the case of two different cosmo- 
logical models = 1 (cases (a) and (b)) and £1 = 0.3 
(cases (c) and (d)). Fig. 3 (a) and (c) show the noise 
free power spectrum (solid lines) the power spectrum mea- 
sured on the reconstructed maps with n = 30 gal/arcmin 2 
(dotted lines) and for fi = 50 gal/arcmin 2 (dashed lines). 
The plateau for the latter two cases is the consequence 
of the intrinsic ellipticities of the galaxies: at small scales, 
the power is dominated by the ellipticity of the galaxies, 
thus P(k) tends to be constant, generally much higher 
than the signal. Figs, (b) and (d) show the difference of 
the power spectrum measurements on the reconstructed 
maps with the noise-free power spectrum; thin dotted 
line is for n = 50 gal/arcmin 2 and thin dashed line for 
fi — 30 gal/arcmin 2 . For visibility, error bars for scales 
larger than 20 arcmin have been dropped. 

In 3D space, these angular scales correspond approxi- 
mately to scales from 1 to 30 h~ 1 Mpc. We leave for later 
studies the problem of inverting the measured projected 
P(k) to the 3D one. This aspect has already been explored 
at large angular scale by Seljak (1997). 

As can be inferred from the power spectra in Figure 
3, the noise has a flat power spectrum, characteristic of a 
white noise process. Fig. (2) illustrates the fact that the 
noise is independent on the underlying k field and follows 
a Gaussian distribution. The noise model introduced by 
Kaiser (1998) based on the weak lensing approximation 
is now compared to the noise found in our simulations. 
The weak lensing approximation applied to Eq. (10) gives 
a local shear estimate 7 = 7 + , where 7 is the true 
shear and = 1/N P e^ s \ the mean intrinsic elliptic- 
ities of N p galaxies in the superpixel p. Since the noise 



This is less evident in the case of non-local mass 
reconstruction. 
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Fig. 1. Example of reconstructions of projected mass maps. The top panels show the initial noise-free k map for either Q = 1 
(left panel) or Q = 0.3 (right panel) with the same underlying linear random field (see Appendix A) and the same rms distortion. 
The bottom panels show the reconstructed k maps with noise included in the shear maps. The maps cover a total area of 25 
degrees 2 . Each pixel has an angular size of 2.5 arcmin 2 and averages the shear signal expected from deep CCD exposures (about 
30 galaxy/arcmin 2 ). The sources are assumed to be all at redshift unity and to have an intrinsic ellipticity distribution given by 
Eq. (44). Such a survey is easily accessible to MEG AC AM at CFHT. The precision with which the images can be reconstructed 
and the striking differences between the two cosmological models demonstrate the great interest such a survey would have. 



components are assumed to be spatially uncorrelated, the 
statistical properties of the noise are, 

(ei s \0 i )$\e j ))=a 2 X f3 5«j, (13) 

where a, (3 — (1,2), S K is the Kronecker symbol and cr £ is 
the variance of one component of the intrinsic ellipticities 
in one supcrpixcl. The shear and the intrinsic elliptici- 
ties of the galaxies are uncorrelated in the weak lensing 
approximation. The measured power (on the noisy mass 
maps) can then be expressed only in terms of the true 
power P K (k) (the one we want to estimate) and the power 
spectrum of the noise (Eq.(13)). If k denotes the Fourier 



transform of the measured convergence, its power spec- 
trum is given by 



(k 2 (k))^nP K (k) + a 2 e . (14) 



This equation is only valid for a compact survey, where n 
is the mean number density of the galaxies per superpixel, 
and a e is the mean value of cr e over the survey. For a sparse 
survey, the first term in Eq. (14) is changed into a con- 
volution term, but the noise contribution to the observed 
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Fig. 2. Histograms of the difference between the noisy mass 
map and the initiai mass map. The vaiues of k have been se- 
lected in bins. The histograms that are found to be all com- 
patible with a Gaussian distribution with a fixed mean and 
width. 



mate of the power spectrum, down to the smallest scales 
considered here. 

The behavior at scales smaller than our pixel size re- 
mains partly an open issue for two reasons: first, due to the 
smaller number of galaxies, the convergence of the recon- 
struction process as well as the stability of the noise prop- 
erties has to be investigated. Second, at small scales the 
gravitational distortion is larger than only a few percent, 
and it can go up to infinity on the critical lines. Therefore, 
estimating the variance of the galaxies intrinsic ellipticity 
distribution, arcs and arclets should be removed. This is- 
sue can only be addressed in high resolution simulations 
like those performed by Jain et al. (1998). 

3.2. Power spectrum cosmic variance 

Although the estimate of the power spectrum described 
above is unbiased, the cosmic variance has also to be ex- 
plored to come up with an optimal observational strategy. 

In estimating the cosmic variance of the power spec- 
trum, Gaussian statistics is usually assumed. This hy- 
pothesis is tested comparing the cosmic variance assuming 
Gaussian statistics to that of the simulated mass maps and 
that reconstructed using the two different level of noise de- 
fined before. 

In the case of Gaussian statistics high order moment 
are related to the second order moments via the following 
relations: 



power spectrum remains independent on that power spec- 
trum. A convenient way to estimate of is to take 



''pix 




(15) 



where n p i x is the number of superpixels and N p the num- 
ber of galaxies in the superpixel p. This is nothing else 
but the variance of the observed ellipticities of the selected 
galaxies 5 . 

The thin dashed and dotted straight lines on Fig. 
3 correspond to the expected noise power spectra (for 
n = 30 gal /arcmin 2 and n — 50 gal/arcmin 2 ) according 
to the Kaiser's model given by Eq. (14). It perfectly fits 
the noise part of the mass reconstructed with noisy data, 
whatever the cosmological model and the noise level. This 
is true even for the open cosmological models for which 
stronger non-linearities could have produced a stronger 
coupling with the noise. The fact that the noise compo- 
nent is pure white noise with an amplitude in agreement 
with the theoretical prediction is a remarkable result since 
the full non-linear equations were used, and it shows that 
the weak lensing approximation can be safely used to re- 
move the noise component and to get an unbiased esti- 



(K(ki)...K(k 2p +i)> = 0, 

v 

(ft(ki).../t(k2 P )) - ]T H(^-i)k(k 2l )), (16) 

perm i— 1 

which means that physically, the frequencies of a Gaus- 
sian field are not coupled, and that the 2p moment at a 
given frequency is only determined by the power at that 
frequency. 

We consider a compact survey of size 0, for which the 
number of modes available at a frequency is maximum. 
Thus, following Fcldman ct al. (1994) and Kaiser (1998), 
the cosmic variance Op (k) of P K (k) is given by the square 
of the measured signal (S 2 (k)) 2 (which depends on <f e ), 
divided by the number of independent modes used to de- 
termine it, AAT(k) in the k-annulus (k, k + Ak). Simple 
modes count gives 



AN(k) = it 



k Ak 
K o 



(17) 



where k = 2tt/Q is the fundamental frequency, thus 
2V^(£ 2 (k)) 



eVFAk 



(18) 



The noise depends only on the variance, not on higher mo- 
ments of the intrinsic ellipticity distribution 



Note that in this hypothesis the cosmic variance in inde- 
pendent on the amplitude of the fluctuations. This is not 
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Fig. 3. Power spectrum analysis of the projected density field. The upper series of plots are for = 1 and the lower series for 
f2 = 0.3. For all the plots, the solid lines show the power spectrum of the noise-free mass maps, before any mass reconstruction. 
The dotted lines correspond to the power spectrum estimation on the reconstructed noisy maps with a number density of galaxies 
of n — 30 gal/arcmin 2 , and the dashed lines for n — 50 gal/arcmin 2 . The left panels show the power spectrum estimates from 
the reconstructed noisy maps, compared to the true power spectrum. The right panels show the power contribution due to 
the noise, it is the difference of power spectra between the reconstructed noisy maps and the true power spectrum. The thin 
dashed and dotted straight lines show the expected value of the noise contribution in the simple linear noise model described 
by Eq.(14). It fits remarkably well the true noise level. 



the case when the non-linear couplings are taken into ac- 
count as it can be seen in Fig. 4. It shows the cosmic vari- 
ance for flat ((a), (c), (e)) and open (((b), (d), (f)) cosmo- 
logical models, (a) and (b) correspond to a 5x5 degree sur- 
vey with z s — 1, (c) and (d) with z s = 1.5, and (e) and (f) 
a 10 x 10 degree survey with z s = 1. On each plot, the thin 
solid line shows the Gaussian cosmic variance, the thick 
solid line shows the true cosmic variance without noise, 
the dotted line the noisy maps with n — 50 gal/arcmin 2 
and the dashed line with n — 30 gal/arcmin 2 . The vertical 
axis gives the error on the power spectrum measured at a 
given scale. 

The departure from Gaussianity appears for scales be- 
low 10', the effect is however more important in the open 
case model for which non linearities are stronger. Open 
models ((b), (d), (f)) give almost the same features as 
for the flat models (((a), (c), (e)), although the cosmic 
variance is smaller. This is clearly a consequence of a 
higher power spectrum signal at low scales for these mod- 
els, which is visible when comparing Fig. 3 (a) (fiat) and 



(c) (open). Thus, as expected, the intrinsic shape of the 
power spectrum affects the cosmic variance (Kaiser 1998). 

Going deep in rcdshift (by comparing (a) and (c), or 
(b) and (d) for the open case) clearly improves the cosmic 
variance at small scale, since the gravitational distortion 
is stronger. However this stronger distortion does not im- 
prove the large scale power estimation because the cosmic 
variance at these scales only depends on the whole volume 
survey. 

If shorter wave vectors are observed, (for a 10 x 10 
degree survey), which is visualized on (e) and (f), a gain 
of 2 is reached at all scales, as a direct consequence of 
global increase of the number of modes in Eq.(17). For 
small scales, from the point of view of the statistics it is 
in fact equivalent to observe deep over a small area than to 
observe shallower over a large area (which reflects the fact 
that the dashed lines in (e) and (f) (n = 30 gal/arcmin 2 ) 
are almost the same as the dotted line in (a) and (b) 
(n = 50 gal/arcmin 2 ). But on the other hand, the shal- 
low large survey gives a better estimate of the power at 
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Fig. 4. Cosmic variance for flat (left panels) and open (right panels) models. The thin solid lines are the cosmic variance 
expected for a Gaussian density field given by Eq. (18). The thick solid lines show the true (measured on the simulations) 
cosmic variance. The dotted and dashed lines show the cosmic variance on the reconstructed mass maps, respectively with 
n = 50 gal/arcmin 2 and n — 30 gal/arcmin 2 . 



large scales than the deep survey. Since these two observa- 
tional strategies require the same total exposure time it is 
clear that wide shallow surveys are better than small deep 
surveys. As it will be shown in Sect. 4, this remains true 
for the high order moments in real space. Moreover, deep 
surveys show more and more distant galaxies for which 
the redshift distribution is more uncertain. 

4. Moments in real space 

4-1. Signature of the normalization and non-Gaussian 
properties 

Weak lensing carries more information than the ampli- 
tude and shape of the dark matter power spectrum. Fig. 1 
demonstrates that distortion maps of the same amplitude 
(and with very similar power spectrum as can be checked 
in Fig. 2) can display very different features. On these 
maps the variance of the local convergence is the same, 
but the amount of non-linearities is very different. For low 
rt universes, the same amount of distortion can be reached 
only with a rather large value of as thus corresponding to 
a much more evolved dynamics. As a result the differ- 
ence between the underdense and the overdense regions 



is more pronounced. The 'voids' tend to occupy a much 
larger area, whereas the super clusters tend to be sharper. 
These features appear because of the non-linear couplings 
contained in the gravitational dynamics. At large scale the 
use of Perturbation Theory has proved to be extremely 
good in predicting the emergence of such properties. All 
these calculations are based on the hypothesis that the 
initial conditions were Gaussian, which we will assume as 
well. 



It has already been stressed (BvWM) that the de- 
parture from a Gaussian statistics is described by the 
skewness of the probability distribution function (PDF) 
of the local convergence. We will restrict our analysis on 
the skewness basically for two reasons: it is beyond the 
scope of this paper to explore all possible indicators of the 
non-Gaussian properties, and we know that the approxi- 
mate dynamics we have adopted reproduces correctly the 
skewness of the local PDF (see Appendix A). In addition 
the lens-lens coupling and the Born approximation terms 
which are known to be small for the third moment are 
probably more important for higher orders, and this re- 
quires a complete dedicated work. 
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Let us summarize the expected results. For a top-hat 
window function we expect to have, 

-(n+2)/2 



a K w 0.01 a 8 Q£ 



(I..S /" 



\ldeg. 



,0.75 



S3 



40 0, 



-0.8 ,-1.35 



(19) 
(20) 



where cr K is the rms value of k at the scale 6q, n is the 
index of the power spectrum, erg is the 3D rms density at 
8/i _1 Mpc scale and z s is the mean redshift of the sources. 
The computed skewness S3 is expected to be independent 
on the normalization of the power spectrum. It is only 
weakly dependent on the shape of the power spectrum as 
well as on the cosmological constant A (see discussion). 
The skewness would then be a very robust way 6 of deter- 
mining the density parameter fl . 

Although the moment analysis is generally performed 
on the basis of a top-hat filter there is a priori no rea- 
son to limit our investigations to this filter. In particular 
SvWJK have proposed the use of an alternative function, 
the compensated filter that might prove more efficient to 
constrain ft, with a lower cosmic variance. 

4-2. Top-hat versus compensated filters 
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Fig. 5. Spectral response for the family of top-hat filters used 
in this work. The finite size effects are visible for the largest 
smoothing scales. The curves are not regularly spaced because 
of pixelisation effects. 



Compensated filters were considered by SvWJK as a 
way to measure the convergence directly from the galaxy 
shape. They use the filter V{0) of size C that defines the 



quantity M ap as 



M ap (0) 



Ja 



V(9) lt {0), 



(21) 



6 It is worth reminding that this is only possible if the redshift 
of the sources are perfectly known (which we assume here) . 
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Fig. 6. Spectral response for the family of compensated filters 
used in this work. 



where j t is the tangential component of the shear field. 
M ap is the convergence field filtered by U(0), 



M ap (0) 



Jo 



d 2 eu(0) k(0), 



(22) 



where U(9) has to be the compensated filter and it is 
related to the arbitrary filter V through 



V{9) 



(23) 



this latter filter is a compensated filter 7 and, for instance, 
a convenient filter to use is given by, 



U{9) 



3 

rife. 



(24) 



The spectral responses (defined as the squared amplitude 
of the Fourier transforms of the window funstions) of the 
filter's family used in this work are shown in Fig. 5 (for the 
top-hat) and by Fig. 6 (for the compensated filter U(9)). 
Clearly a field smoothed with a top-hat filter of size C 
is sensitive to fluctuations of size larger than C that con- 
tribute also to the cosmic variance of the moments. On the 
other hand, a compensated filter integrates the fluctuation 
modes only around the target frequency C and any power 
at lower or larger scales will affect neither the signal nor its 
cosmic variance. SvWJK showed that the de-correlation 
properties of compensated filters are by far better than 
for top-hat filters because any power at small wavelengths 
between two disconnected fields is highly suppressed (see 
Fig. 8 of their paper). Thus the cosmic variance should be 
smaller for a compensated filter than for a top-hat filter. 

7 A compensated filter is a filter with a vanishing mean value. 
This explains why M ap can be measured only from the distor- 
tion field only since the mass-sheet degeneracy is removed by 
the used type of filters. 
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However this attractive feature comes with a price: a com- 
pensated filter needs to be sampled by a larger number of 
galaxies. In other words, the shot noise has a larger effect 
on the aperture mass M ap than on the top-hat filtered 
mass at the same scale. A compromise has to be found, 
that depends on the functional shape of the compensated 
filter and on the shape of the power spectrum. 

4-3. Moment estimations and shot noise corrections 

Due to the intrinsic ellipticities of the galaxies and the cos- 
mic variance (see for instance Szapudi & Colombi 1996, 
Colombi et al. 1998), estimates of the moments of k from 
the reconstructed map are biased. We show here that the 
shot noise can be accurately calculated and the estimated 
moments corrected. It is worth noting that SvWJK has 
shown that we can find an unbiased estimator of the mo- 
ments of M ap , which completely cancel the shot noise cor- 
rection problem. It is based on the measurement of the 
tangential shear 74. Unfortunately, the measurable quan- 
tity is gt (the tangential reduced shear) rather than j t , 
and unless this is taken into account, the estimator given 
in SvWJK is no longer unbiased. In other words, the shot 
noise correction problem is shifted to an estimator correc- 
tion problem. 

It was shown in Section 3.1 that the shot noise leads to 
a pure white noise in the reconstructed convergence maps. 
The amplitude of this noise can be obtained by measur- 
ing the observed ellipticities of galaxies as described for 
the power spectrum estimation. Therefore, estimates of 
the variance and the skewness of the convergence in (19) 
and (20), corrected from the intrinsic ellipticities of the 
galaxies, are obtained by simply removing the noise term 

((o s ') ) in the second moment. Note that the skewness 
correction only requires the correction of the variance since 
the third moment is not affected by the noise. As the ana- 
lytical calculation of the noise term for compensated filters 
can be rather cumbersome, it is estimated using Monte- 
Carlo simulations. 

Even after that noise correction, finite sample effects 
may bias the estimations of the second moment and of S3 
and may increase significantly the cosmic variance. This 
difficulty was partly investigated in BvWM with the use of 
perturbation theory. They pointed out that the accessible 
geometrical averages are expected to be smaller than the 
true ensemble averages, and that a dispersion is expected 
in the measurements (cosmic variance): 

— the bias that affects the expectation values was found 
to be proportional to the variance at the sample size 
divided by the one of at the filtering scale; 

— the scatter was found to be proportional to the rms of 
k at the sample size. 

These estimates were done fully in perturbation the- 
ory, with numerous approximations (in particular it was 
assumed that the sample size was much bigger than the 



smoothing scale which is probably an erroneous approxi- 
mation for most of the cases considered here). For accurate 
investigations of all these effects that take into account 
both the Poisson noise and the finite volume effects see 
Szapudi & Colombi (1996). 

4-4- Results 

We now turn to the measurement of moments in the simu- 
lated fields. The same simulations used for the power spec- 
trum analysis are used here. The results are given in Figs. 

7 (variance (k 2 ) in a 5 x 5 degree field with z s = 1) and 8 
(skewness s 3 in a 5 x 5 degree field with z s = 1). For each 
of these figures the plots are organized in the same way: 
the first raw ((a), (b) and (c) plots) is for the flat model, 
the second raw ((d), (e) and (f) plots) for the open model. 
It corresponds to the first two raws of Table 1. The first 
columns (a) and (d) show the estimator measured with a 
top-hat filter, the second columns (b) and (e) with a com- 
pensated filter, and the third columns (c) and (f) show the 
signal to noise ratio of these estimators. In plots (a) , (b) , 
(c) and (d) the solid lines give the estimators measured in 
the noise-free k maps, the dotted lines in the noisy recon- 
structed maps with n = 50 gal/arcmin 2 and the dashed 
lines with h = 30 gal/arcmin 2 . The dotted-dashed lines 
show the estimators measured on the reconstruction with 
noise with n — 30 gal/arcmin 2 corrected from the noise. 
Since the case n = 50 gal/arcmin 2 gives the same results 
they are not plotted. On the signal to noise plots (c) and 
(f) the thin solid lines show the results for a top-hat filter 
and the thick solid lines for a compensated filter. The re- 
sults obtained for the noisy maps with n = 30 gal/arcmin 2 
or n = 50 gal/arcmin 2 corrected from the noise are respec- 
tively given by the dashed and dotted lines (either thin or 
thick). 

4.4.1. Noise correction 

The noise correction as described in Sect 4.3 gives un- 
biased results, as it is expected for superimposed white 
noise. This is true even when the correction is two orders 
of magnitude higher than the signal (see for example Figs. 

8 plots (b) and (e)). This confirms the simple properties 
of the noise in the reconstructed mass maps already found 
in Sect. 3.1. 

4.4.2. The variance 

The variance obtained for a top-hat filter is slowly decreas- 
ing with an increasing scatter, as expected. For the com- 
pensated filter the curves are almost flat as expected from 
the shape of the power spectrum. The noisy maps (dotted 
and dashed lines in Fig. 7. (a) and (b) display higher val- 
ues for the variance. Once it is corrected, the results are in 
perfect agreement with the noise- free simulations. The sig- 
nal to noise ratio is basically not affected by the shot noise 
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Fig. 7. The measured variance of the convergence, and the corresponding signal to noise ratio. The upper panels show the fi = 1 
case, while the bottom panels correspond to the £1 = 0.3 case. The left and middle panels are respectively the measured variance 
with a top-hat and a compensated filter. On these plots, the thick solid line is the true variance measured on the noise-free 
maps, the dashed line and the dotted line for the noisy reconstructed mass maps (with respectively n = 30 gal/arcmin 2 and 
n — 50 gal/arcmin 2 ). The dotted-dashed line is the variance measured from the n — 30 gal/arcmin 2 case and corrected from the 
noise. The right panels show the signal to noise ratio of the variance detection with the top-hat (thin solid line) and compensated 
(thick solid line) niters. Dotted and dashed lines have the same meaning as for (a), (b), (d) and (e) but here the variance has 
been corrected from the noise. 



for a top-hat filter, it shows that going deep does not im- 
prove the measurement precision. The compensated filter 
reveals much more sensitive to the shot noise as predicted 
in section 4.2 since we can see on plot Fig. 7. (c) and (f) 
(thick lines) the bell shape of the signal to noise ratio, 
with a significant reduction of the measurement precision 
at the smallest available scales. The open and flat cases 
show basically no differences. It can be seen however that 
the signal to noise ratio for a compensated filter is sightly 
lower for the open case. We interpret this effect as due 
to the presence of more nonlinear couplings in the maps. 
Remarkably, the precision with which the variance can be 
measured in some specific k range reaches 5%. 



4.4.3. The skewness 



The skewness can be accurately measured at the smallest 
scales (in Fig. 8 only the error bars for the first four points 
have been drawn). The skewness is decreasing with scale in 
the two cases. Once again, the noise correction applied to 
the reconstructed maps allow to recover the skewness with 
a surprising accuracy. The signal to noise of the skewness 
is still not affected by the noise for a top-hat filter, while 
for the compensated filter the situation is worse than for 
the variance; for instance, at the smallest scale, the signal 
to noise is almost one order of magnitude smaller on the 
reconstruction with noise than on the noise-free maps for 
the = 1 case. The two models, open and flat, provide 
us with very different magnitude for S3. Their skewness 
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Fig. 8. Same as Fig. 7, but for the skewness of the convergence, S3. 



ratio is 3 as expected from perturbation theory 8 , with a 
significance of the separation at roughly 6o\ This is ob- 
served in case of the top-hat as well as the compensated 
filter and this confirms the fact that the skewness of the 
convergence can strongly separate low and high density 
universes. 

To be more precise we present the actual histograms 
of the measured skewness in Fig 9 which demonstrates 
clearly that the two cosmologies can be easily separated. 
One can see that the scatter in S3 is roughly the same in 
the two cases and that the difference in the relative pre- 
cision is due to the differences in the expectation values. 
This plot also shows that the distribution of the measured 
S3 is quite Gaussian. 
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To reduce the cosmic variance, we reject the pixel values 
that are above 4 a. It slightly reduces the value of S3 for the 
open case. 



Fig. 9. Histograms of the values of S3, top- hat filter, for $1 = 1 
(solid lines) and Q = 0.3 (dashed lines) for a 5 x 5 degree 
survey (thick lines) and a 10 x 10 degree survey (thin lines). 
The angular scale is the pixel size 2.5'. 
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4-5. Comparisons of different observational scenarios 

The results presented previously had been obtained from 
the full mass reconstruction of 240 maps 9 as described 
in Appendix B. Since it was demonstrated in the preced- 
ing section that noise acts as a pure de-correlated white 
noise in the reconstructed n maps, we pursue our analy- 
ses of large series of simulated helds simply by adding the 
noise on the initial n maps (especially for 10 x 10 degree 
data sets for which convergence reconstruction would take 
typically one and a half hour on DEC PWS-500 comput- 
ers). The subsequent analysis are therefore made with this 
simplified scheme. 

To complement the previous cases, we have built and 
analyzed the cosmic variance on 60 maps for each of the 
following models: open (SI — 0.3) and flat (Cl = 1) cos- 
mologies, a survey size of 5 x 5 degrees for z s = 1.5, a 
survey size of 10 x 10 degrees for z s = 1, with the power 
spectrum of Eq. (38) and a survey size of 5 x 5 for z s = 1 
with a CDM spectrum. 

4.5.1. Effect of the survey size 

By increasing the total area of a factor 4 we increase the 
signal to noise on the variance and the skewness with a 
top- hat filter by almost a factor 1.7. This can be seen 
in Fig. 9 when comparing the 5x5 degree case with the 
10 x 10 case. With a compensated filter, the signal to noise 
ratios of the variance and of the skewness are increased by 
exactly a factor 2, thus improving more rapidly with the 
sample scale than the top-hat window function. This is 
expected from the de-correlation properties of those filters. 
It makes the compensated filter actually more attractive 
for such a large survey. 

4.5.2. Effect of the source redshift 

Fig. 10 shows the effect of a change in the mean source red- 
shift. For the galaxies that are further away, the variance 
of k is larger since the gravitational distortion is stronger, 
conversely the skewness is smaller since the accumulated 
material along the line of sight creates a field that is more 
and more Gaussian. The surprising result is that the sig- 
nal to noise of these quantities does not depends strongly 
on the redshift of the sources. This means that it will be 
a waste of time to observe at high redshift, while it will 
basically not improve the precision of the measurement. 
Things are slightly improved for the compensated filter, 
but fundamentally, the results with the top-hat filter show 
that we do not learn more by increasing the redshift. Note 
that there is no improvement due to the increasing galaxy 
number density if the survey size is unchanged (see section 
4.3). In addition, high redshift surveys may create new 

9 This corresponds to the two noise levels for each cosmolog- 
ical model, with 60 mass maps for each case. 



problems such as uncertainties due to the Born approxi- 
mation, the lens-lens coupling or the recently investigated 
source clustering effect (Bernardeau 1998). 



4-6. BG versus CDM power spectrum 

In order to test the robustness of the skewness as an esti- 
mator of Cl independent on the power spectrum we re-ran 
our simulations for a standard CDM power spectrum. Fig. 
11 shows the comparison between the CDM model (dot- 
ted line) and BG power spectra, which clearly shows that 
CDM contains more structures at small scale by looking 
at the variance plots (a), (b). As predicted in BvWM, 
the skewness of the convergence if almost unaffected by 
the change of power spectrum (see Fig. 11 (d), (e)), but 
there is a small improvement in the signal to noise for the 
flat model (because of the larger power at small scale for 
CDM). On the other hand, the variance is strongly af- 
fected and in particular there is a significant decrease of 
power at large scale compare to BG spectrum. Note that 
the compensated filter yields a more accurate representa- 
tion of the underlying power spectrum than the top-hat 
filter (because it is a pass band filter), and leaves the an- 
gular dependence of the variance unchanged compared to 
the spectrum Eq.(38) excepts at large scales. 

The skewness is thus a robust estimator of Cl fairly 
insensitive to the power spectrum. 



0.12 T 
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Fig. 12. Histograms of the values of S3 for 5x5 degree survey, 
top-hat filter, for Cl — 1 for a BG spectrum with as = 0.6 and 
z s = 1 (thick solid line), as = 0.6 and z s = 1.5 (thin solid line), 
a s — 1.0 and z s = 1 (dotted line) and a CDM spectrum with 
as = 0.6 and z s = 1 (dashed line). The angular scale is the 
pixel size, 2.5'. 
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Fig. 10. For a flat cosmological model, comparison of different observational strategies between a survey of size 10 x 10 degrees 
(dashed lines) a 5 x 5 degrees survey for sources at mean redshift 1.5 (dotted lines) and a 5 x 5 degrees survey for sources at 
mean redshift 1 (solid lines). Left panels are for the top-hat filter and middle panels for the compensated filter. The thicker 
lines in the right panels hold for the compensated filter. 



4-7. Effect of the normalization 

The skewness S3 is found to be independent on the normal- 
ization, as expected from the Perturbation Theory. The 
signal to noise ratio however is increased by about 40% in 
the case of high normalization ag = 1. These results are 
summarized in Fig. 12 that shows the histograms for vari- 
ous cosmological cases. It demonstrates that the skewness 
is clearly independent on the shape and normalization of 
the power spectrum. However there is a strong dependence 
on the mean redshift of the sources. If the signal to noise 
ratio for s 3 depends on the cosmology it is independent 
on the mean source redshift. This again is favoring rather 
shallow surveys. 

4-8. Beyond the skewness to measure SI? 

The skewness of the PDF of the local convergence does 
not entirely characterize the PDF itself. Thus it is natural 
to measure higher order moments of the convergence to 
probe the cosmology. 



It is clear that the skewness breaks the degeneracy be- 
tween the power spectrum and the cosmological param- 
eters, and is completely insensitive to the normalization. 
Bcrnardeau (1995) already noticed that in the case of cos- 
mic density field or cosmic velocity field, the ratio Si/s\ 
calculated from the perturbation theory with a top-hat 
filter is almost a constant, and independent on the un- 
derlying cosmological model. This work was recently ex- 
tended to the lensing case (Bernardeau, 1998) where he 
found s 4 (k)/s 3 (k) 2 ~ 2. If all systematics of the gravita- 
tional lensing measurement can be controlled, search for 
such a magic number in our Universe would be a strong 
indication of validity of the paradigm of the gravitational 
instability scenario started from Gaussian initial condi- 
tions. On the other hand S4 may be a new way to measure 
the density parameter (but not totally independent on the 
skewness). In our maps we find that the noise correction 
still works for the kurtosis, and that a compensated filter 
is more efficient than the top-hat filter (at least for noise- 
free data). In addition the kurtosis appears to be a fairly 
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Fig. 11. Comparison of different choices of power spectra for a 5 x 5 degree survey and a flat cosmological model. The dashed 
lines correspond to an SI = 1 CDM spectrum with a s = 0.6 and z s = 1, dotted lines to a BG spectrum with as = 1 and z s — 1, 
and the solid lines to a BG spectrum with as = 0.6 and z s = 1. 



good discriminant for Q. Unfortunately the signal to noise 
ratio remains lower than for the skewness, which make it 
more difficult to measure. Moreover, it is more sensitive 
to the usual lensing approximations (Born approximation 
and lens-lens coupling terms) as well as source clustering 
(Bcrnardeau 1998). The error bars of the kurtosis found 
in our simulations are so large with a 5 x 5 degree survey 
that it is impossible to measure it at scales larger than a 
few arcminutes. It turns out that the ratio S4(k)/s3(k) 2 
is ~ 1.4 for the flat case and ~ 1.8 for the open case 
for a top-hat filter, while it is ~ 5.6 for the flat case and 
~ 3.4 for the open case for a compensated filter. It is not 
the scope of this work to compare further the differences 
between the two filters, but we want to point out that 
filtering may be an interesting way to change the depen- 
dence of an estimator versus the cosmology. This should 
be studied in order to search for optimal measurements 
of higher order moments. More generally, the whole shape 
of the PDF could probably be used with more efficiency 
than the skewness alone. In a regime of small departure 
from a Gaussian distribution it can for instance be fruit- 



ful to describe the shape of the PDF with an Edgeworth 
expansion that takes into account the first few moments 
(Juszkiewicz et al. 1995, Bernardeau & Kofman 1995). For 
instance, from the Edgeworth expansion it is easy to show 
that the fraction of values of the convergence that is above 
the average value, P(k > ( < k) >), is 

"^^-i < 25 » 

The results obtained from this formulae are in good agree- 
ment with those obtained from the direct measurements, 
however with a slightly larger cosmic variance. 

Finally, the non-Gaussian features can also be char- 
acterized with topological indicators (which have been 
shown to be fruitful for the analysis if CMB data already, 
see for instance Winitzki & Kosowsky 1998, Schmalzing & 
Gorski 1997). What cosmic variance could be derived from 
the joint use of topological quantities or/and information 
on the shape of the PDF is left for further investigations. 
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5. Discussion 

In this paper, we have focussed our investigations on scales 
larger than 2.5' thus allowing ourself the use of a simpli- 
fied dynamics. This is a complementary approach to the 
ongoing investigations by Jain et al. (1998) who analyzed 
the high order moments and the power spectrum of the 
calculated from ray tracing in high resolution simulations. 
Although they only analyze the statistical properties of 
the convergence without noise, their approach completes 
our own towards the smaller angular scales. They have 
shown in particular that the skewness of the convergence 
is significantly higher at small scales (0.1 arcmin) than 
the theoretical expectations of perturbation theory due to 
highly non-linear structures, (this was already mentioned 
by Gaztanaga & Bernardeau 1998 and whether this be- 
havior can still be described by means of perturbation 
theory with the help of loop correction terms is still an 
open question). Unfortunately, even at these small scales, 
they are not able to analyze the cosmic variance because 
of the small number of realizations. The use of high res- 
olution simulations probably prevents a detailed analysis 
of this quantity. 

The preceding sections provide quantitative estimates 
of the capability of mass reconstructions from weak lens- 
ing measurements at large scale to probe the large scale 
structures, as well as the cosmological parameters. We 
have shown how the projected mass distribution can be 
reconstructed accurately from the observed shape of the 
galaxies. Two complementary analysis have been exam- 
ined, the power spectrum and the non-Gaussian features 
through the high order moments. 

For the power spectrum estimation, the best sampling 
strategy (i.e. the question of sparse or compact surveys) 
is not discussed in this paper (see Kaiser 1998 for a dis- 
cussion), but our results show that in order to probe the 
smallest scales of mass fluctuations a deep (but narrow) 
survey is required that diminishes the cosmic variance 
caused by the shot noise. Once this is done, the survey 
can be extended in a more shallower manner to probe 
the power spectrum at scale where the cosmic variance 
caused by shot noise becomes unimportant. At this stage, 
the question of a sparse or compact survey is a matter of 
choice, depending on scientific interests. 

Concerning the moments measurement, a summary is 
given in Table 2 (for a top-hat filter) and 3 (for a com- 
pensated filter) which shows the smallest accessible error 
on the measurement of the variance and the skewness of 
the convergence for different observational contexts. This 
demonstrates that weak lensing measurements can reach a 
few percent precision on fl for a reasonable survey size. We 
pointed out that the results are not much deteriorated for 
a number density of galaxies of 30 gal/arcmin 2 compared 
to 50 gal/arcmin 2 , whereas the realization of the survey 
in the latter case requires a factor 3 more observing time. 
It suggests that large shallow surveys would be more ade- 



Table 2. Relative error on the measured variance and skewness 
of the convergence in different observational contexts (in tables 
ni = 30 gal/arcmin 2 and 712 = 50 gal/arcmin 2 ) for a top- hat 
filter. The numbers correspond to the smallest error within the 
range of scales considered in this work. 



Observational 
constraints 


Q = 1 


Q = 0.3 


Scr K /a K 


5S 3 /S 3 


S<j k /<j k 




5x5 deg 
z s = 1 


m 


0.051 


0.190 


0.047 


0.118 


ri2 


0.050 


0.186 


0.047 


0.105 


5x5 deg 
z s = 1.5 


m 


0.043 


0.246 


0.038 


0.127 


712 


0.043 


0.223 


0.038 


0.123 


10 x 10 deg 

Z s — 1 


m 


0.029 


0.116 


0.024 


0.060 


712 


0.029 


0.110 


0.023 


0.054 



Table 3. Same as Table 2 for a compensated filter 



Observational 
constraints 


fi = 1 


n = 0.3 


8a K 1 a K 


5S 3 /S 3 


5<j k / ' o K 


SS3/S3 


5x5 deg 

Z a = 1 


711 


0.028 


0.391 


0.044 


0.209 




0.021 


0.310 


0.043 


0.174 


5x5 deg 
z s = 1.5 


m 


0.016 


0.341 


0.030 


0.193 


712 


0.016 


0.310 


0.029 


0.172 


10 x 10 deg 

Z 3 — 1 


m 


0.014 


0.171 


0.020 


0.110 


712 


0.014 


0.140 


0.020 


0.074 



quate since it efficiently reduces the cosmic variance. Such 
a strategy would be much more comfortable with respect 
to some systematics like the redshift distribution of the 
sources (which can be determined easily if the sources are 
closer) , the source clustering effects that are limited if the 
source distribution is narrow (Bernardeau 1998), and pos- 
sibly the morphological evolution of distant galaxies, in 
particular if most distant galaxies are composed of many 
merging substructures. 

Note that there is still room for potential improvement 
of the signal to noise ratio we are obtaining. Indeed, since 
we are limited by construction to simulations of scales 
larger than 2.5' we do not know whether the cosmic vari- 
ance can be reduced by observing the moments at smaller 
scales, as it is suggested by the results obtained with a 
top-hat filter (for which the signal to noise curves never 
bend down at small scales) . The optimal size for the mea- 
surement of the variance and the skewness might in fact 
correspond to the arcmin scale. However, there are sev- 
eral issues that are generally believed to be irrelevant for 
weak lensing at large scales, but are probably potential 
difficulties at small scales: 

— Our work implicitly assumes a constant power spec- 
trum below the pixel size. A mass reconstruction from 
real data should include small scale features such as 
cluster lensing, and the propagation of the noise from 
this peaks of signal is not known. Moreover that abil- 
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ity of the x 2 method to match the noise properties 
correctly at small scales should be reinvestigated. 

— Born approximation and lens-lens coupling terms are 
stronger at small scales (see SvWJK). A quantifica- 
tion of this effect on simulations would be necessary to 
decide at which scale the measurement of moments is 
optimal. 

— Source clustering (investigated at large scales by means 
of perturbation theory by Bernardeau 1998) could also 
have a significant impact at small scale. 

The possible systematics caused by all these effects 
have to be investigated in order to estimate the precision 
of weak lensing surveys. 

Moreover the moments are not necessary the best 
means to distinguish between different cosmological mod- 
els (not to mention the fact that the results are anyway 
sensitive to the choice of the filter). For example, it seems 
that by comparing open and flat universes on Fig. 1 topo- 
logical tools should be as well a strong discriminant of 
cosmological models. The statistical instruments are ex- 
tremely diverse. Among the possible tools the results of 
multi-scale filtering, as provided by wavelet transforms, is 
likely to be a good candidate. 

6. Conclusion 

We have investigated in detail how weak lensing obser- 
vations could be used to measure the projected power 
spectrum and to discriminate among different cosmo- 
logical models. In order to analyze the cosmic variance 
of the results (which requires a large number of real- 
izations), we used a two dimensional, second order La- 
grangian dynamics to generate the convergence fields. 
In each of the observational situation and cosmological 
model, 60 different maps has been generated. Cosmologi- 
cal models include: flat model (fi = 1), open model (CI = 
0.3), Baugh & Gaztahaga initial power spectrum, stan- 
dard CDM spectrum. Observational contexts include: high 
noise (30 gal/arcmin 2 ) and low noise (50 gal/arcmin 2 ) 
level, "small" survey size (5x5 degrees) and large sur- 
vey size (10 x 10 degrees), two different mean redshift of 
the sources (z s = 1 and z s = 1.5), and the use of top-hat 
and compensated filters. From these convergence maps, 
a shear map is calculated, on which the noise is intro- 
duced according the selected observational context, the 
mass map is finally reconstructed (using a \ 2 method) 
and analyzed. 

Our results are: 

— At scales larger than 2.5' the \ 2 reconstruction method 
is a very stable process which does not produce any 
boundary effects, spurious signal, and which leaves the 
noise properties unchanged (the noise does not prop- 
agates and remains close to the theoretical prediction 
using the linearized lens equation). This permits to 
work with the convergence (which is the physical field 



of interest) instead of the shear, with no loss of infor- 
mation. 

— The shot noise contribution can be removed simply by 
measuring the observed ellipticities of the galaxies, and 
this leads to unbiased estimates of the power spectrum 
and the moments of the convergence (at least up to the 
kurtosis). 

— The precision obtained on the normalization can be as 
low as 2% with survey of 10 x 10, and 5% for f2 (see 
tables 1 and 2) for the most favorable cases (i.e. low 
Vl). 

— The larger shallower surveys are more promising to 
recover the cosmological quantities. In addition to 
that, the redshift distribution of the brighter galaxies 
is known better and their light distribution is more 
regular, making ellipticity measurement from their 
weighted second moment easier and more relevant. 

— The compensated filter yields by far smaller cosmic 
variance than the top-hat filter for the variance of the 
convergence, while it gives unsatisfactory results for 
the skewness. However the cosmic variance decreases 
more rapidly with an increasing survey size for a com- 
pensated filter. 

The MEGACAM project offers the possibility to per- 
form these observations. Indeed the large field of the CCD 
device (1 square degree) and the high image quality at 
CFHT (extended to the edges of the field with the future 
field corrector) provides the ideal instrument to perform 
such a scientific program. 

In the context of the rapid changes occurring in obser- 
vational cosmology it is worth stressing that weak lensing 
surveys offer precious complementary views of our Uni- 
verse and unique tools to probe directly the dark matter 
and to compare with the light distribution at any scales. 
The perspective of determining the projected power spec- 
trum independently of biases is indeed attractive. More- 
over the possibility of determining Qq in a way which is 
independent on the power spectrum, and independent on 
all the methods that have been suggested so far, is also 
extremely precious. We remind that this determination 
relies only on dynamical effects assuming that the large- 
scale structures originate from Gaussian initial conditions 
through gravitational instabilities. In Fig. 13 examples of 
constraints in the (£1, A) plane with weak lensing as de- 
scribes in this work (shaded areas) are presented 10 to- 
gether with the location of the major constraints that are 
expected to be brought either with CDM experiments, 
(constant curvature density, solid lines), or from Type- 
la supernovae measurements (the other straight lines, de- 
scribing constant go values). 

Future CMB experiments can determine the cosmolog- 
ical parameters with a remarkable accuracy, but only when 
some prior is put on the shape of the initial power spec- 

10 The A dependence is taken from the theoretical results given 
by BvWM. 
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0.0 0.5 1.0 1.5 




0.0 0.5 1.0 1.5 



Fig. 13. Constraints that can be brought by weak iensing sur- 
vey in an fl — A plane. The grey areas are the location of 
the 1 and 2-cr bands (respectively darker and lighter bands) 
allowed by a measured skewness that would be obtained with 
either Slo = 0.3 (left bands) or Qo — 1 (right bands). The 
solid straight lines corresponds to a zero curvature universe, 
and the dot-dashed lines to a fixed acceleration parameter, qo- 
The panels correspond to survey of either 5x5 (top) or 10 x 10 
degrees (bottom). 



trum. The curvature will probably be determined with a 
good accuracy in the near future from the position of the 
first Doppler peak on the C; curves. But it appears that it 
can be very difficult to disentangle flo from A (see for in- 
stance Zaldariagga, Spergel & Seljak 1997) from the mere 
temperature (and polarization) fluctuations. This might 
be only possible from a very detailed analyses and the 
power spectra which require not only a good understand- 
ing of the possible systematics that may affect the mea- 
surements, but also specific hypothesis on the regularity of 
the primordial power spectrum. The requirements to mea- 
sure Qq from weak Iensing survey are less strict provided 
the instrumental systematics can be controlled correctly. 

Moreover weak Iensing surveys are able to constraint 
A if it is possible to select efficiently several populations 
of sources (see Villumsen 1996, BvWM for clues of such 
possibilities). We thus think that weak Iensing survey will 
be a major mean for probing the global cosmological pa- 
rameters. It is complementary to the CMB experiments 
and indispensable for breaking the parameter degeneracy. 

The feasibility of the weak Iensing by large scale struc- 
ture program and The full scientific exploitation of the 
weak Iensing effect is reduced to the capability to control 
instrumental systematics, where there are still open cru- 
cial points. The most dedicated problem comes probably 
from the need of a PSF correction able to avoid artificial 
large scale coherent alignments of galaxies. Related to this 
problem is the pixelisation effect. Galaxies shapes are in- 
deed determined from a finite number of pixels with the 
possible introduction of errors and biased estimation on 
the shapes. 

These effects can all be investigated independently. 
Taking several images of the same portion of the sky with 
a shifted position for the camera, or with different cameras 
would definitely test the robustness of the observed dis- 
tortion maps. Moreover we have numerous statistical tests 
at our disposal that can be done on the maps: comparison 
of the 2-point correlation function of the distortion with 
the one of the shear field or with the one of the K field, 
or even by the investigation of quantities that are a priori 
sensitive to the systematics in a different manner, such as 
the correlation of the orientations of the distortion field. 
That all these quantities have related properties is some- 
how due to the fact that in the thin lens approximation 
the 7 field should be curl-free (see for example Luppino 
& Kaiser 1997 for a used of these tests). This is true only 
when lens-lens couplings are neglected but it should be 
possible to investigate those effects in numerical simula- 
tions. 
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Appendix A: Second order Lagrangian dynamics 

In order to investigate the statistical properties of the re- 
constructed mass field that contains a significant amount 
of non-linearities, a non-linear evolution model of large- 
scale structures is required. Second order Lagrangian per- 
turbation theory has the advantage to be extremely fast to 
compute, and rather accurate compared to N-body simu- 
lations. It gives correct values for the skewness 5*3 and we 
expect it to gives correct estimation of the cosmic variance 
(see Munshi et al. 1994, Bernardeau et al. 1994, Bouchet 
et al. 1995). 

Actually we do not even need to do 3D simulations. At 
the level of perturbation theory it is equivalent to perform 
2D second order Lagrangian evolutions of the structures 
on an initial linear map of the projected mass fluctuations. 

Construction of the initial linear map 
The local convergence map is given by 

«(V) = j d Xs n(xs) j dx 

Vk(Xs-x) V k ( X ) 5(V k ( X )<P,x) 



Vk{Xs) 



(26) 



where Xh is the horizon distance, fl the density parame- 
ter, a the expansion factor and n{Xs) the number density 
of sources as a function of the distance x- F° r clarity we 
introduce the lens efficiency function 

w(x) = ~ / d Xs n(x s ) , . , (27) 

so that the local convergence is simply given by the inte- 
gral over the line of sight of, 



K(ip) 



/•xh 
1 ' 



dx w(x) 5(V K ( X )tp,x)- 



The projected density contrast would be, 



= — with ui 



-J 

Jo 



dX w(x)- 



(28) 



(29) 



It is important to keep in mind that the local convergence 
is related to the actual density contrast with a constant 
that depends on the cosmological parameters. In the fol- 
lowing we will assume that all sources are at the same 
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rcdshift (it is not realistic but of no consequences on our 
results here). 

In the linear approximation the k field is expected to be 
a 2D Gaussian field, characterized by a power spectrum, 
P K . with (see Kaiser 1992, 1996, and SvWJK) 11 , 



XH 



dx w(x) 2 P& 



,X 



(30) 



as a result of the relation between k and the 3D density 
contrasts. The initial conditions for k are therefore gen- 
erated by a 2D Gaussian field in Fourier space where the 
complex random variables verify 



(K(k)K*(k')> 



(2n) 2 S D (k- 



k')^(k), 



(31) 



The 2D second order Lagrangian dynamics 



In order to introduce a significant amount of nonlineari- 
ties in the maps, we apply the 2D second order Lagrangian 
dynamics to the projected density, &id- Following the no- 
tations of Bouchet et al. (1992), transposed in the 2D case, 
let us write the 2D Eulerian coordinates 0(q) as a pertur- 
bation series over the 2D displacement field D, 



0(q) = q + eD«(q) + e 2 D< 2 >(q) + o(e 3 ), 



(32) 



where q is the angular Lagrangian coordinate and e a small 
dimension-less parameter. The first order term reduces to 
the Zel'dovich (Zel'dovich 1970) approximation. The di- 
vergence of the second order displacement field can be 
written in terms of the first order solutions, 



Dg(q) = | 



D (1) D (1) 



(33) 



This result is exact for an Einstein-de Sitter Universe. The 
values of the coefficient 3/7 is only slightly altered (about 
1%) for other cosmological models (Bouchet et al. 1992, 
Bernardeau 1994) and in the following we did not take 
into account this dependence. Once the second order dis- 
placement field has been computed, the local 2D density 
contrast can then be written in terms of the Jacobian of 
the transform between the Lagrangian coordinates and the 
Eulerian coordinates, 



<M0(q)) = 



1 



■/(q) 



(34) 



where J 

by 



\d9i/dqj\. The local convergence is then given 



K(0(q))=57 522>(0(q)). 



(35) 



The linear density maps are built on a regular grid from 
random modes following a given power spectrum. The dif- 
ferent quantities up to the Jacobian are generated on this 



Taking H = c = 1. 



same grid by successive Fast Fourier Transforms. There 
is then a technical difficulty to solve in order to get the 
resulting values of k on a regular grid. This is done via 
a local triangulation and an interpolation of the values 
(from standard IDL packages). Note that the continuity 
equation provides us with k and not with the projected 
potenticl. The latter will have to be computed from a sub- 
sequent Fourier transform (see Appendix B). The ampli- 
tude of the fluctuations is such that the displacement field 
does not induce shell crossings 12 . Finally, bands of suffi- 
cient width along the edges were cut out to avoid edge 
effects induced by the displacement filed. 

The skewness with this approximate dynamics 

The skewness of the projected density is given by the skew- 
ness of the 2D dynamics (Munshi et al. 1997), that is, 



density _ 36 3 dlogcr 2 (6>) 



(36) 



J 7 2 dlog<9 

The skewness for the convergence would then simply be 
1 /36 3dlogcr 2 (6») 



convergence 



uJ V 7 2 dlog<9 



(37) 



This result has to be compared with the results obtained 
in BvWM. Their general formula (67) contains some extra 
geometrical factors of order unity (coming from a differ- 
ent averaging procedure along the line of sight). Eq. (37) 
actually corresponds to the approximate form of Eq. (75) 
of BvWM. 

Shapes and normalizations of the power spectra 

The 3D power spectrum given by Baugh & Gaztanaga 
(1996) (BG spectrum) is used in most of our simulations, 



P(k) oc 



[l + (fc/fc c ) 2 ] 3/2 ' 



(38) 



where k c — 0.05 ft-iooMpc - . In a series of simulations we 
compared this model also with the standard CDM spec- 
trum. 

In most cases the fluctuations are normalized accord- 
ing to the convergence field a K = Cste. The value of as 
that has been chosen thus depends on Q in such a way 
that it is equal to 0.6 for a flat Universe (following the 
normalization inferred from galaxy cluster counts (Eke et 
al. 1996, Oukbir & Blanchard 1997). As a result we take, 



0.6 



g(go = 1) 
w(fio) 



(39) 



However, in the open cases it very rarely happened that k 
reached unrealistic large values. The amplitude of k has there- 
fore been arbitrarily limited to 0.1. At most a handful of pixels 
have been affected by this cut-off. 
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It implies that the value of as grows for low values of Q. 
Note that since ZJ ss fi -0 8 , this growth is only slightly 
more important in our case than for the galaxy cluster 
counts (the exponent would be about 0.5 to 0.6). 

Appendix B 

This appendix gives the details of the reconstruction algo- 
rithm, how the g map is generated from an initial n map, 
and the noise model. 

Reconstruction algorithm 

The lensing quantities of interest are defined in the 
Eqs.(6), and the observable is the reduced shear g = 
7/(1 — k). The reconstruction problem is how to infer the 
k map from an observed ellipticity field e b s , knowing that 
(e 6s) = g is an unbiased estimate of the reduced shear? 
Note that this problem is under-constrained since a change 
in the potential tp of the form, 



tp -> Xtp + 



(1-A) 



(40) 



where A is a constant, leaves the reduced shear g invariant, 
but will transform the convergence as k — > \k + 1 — A 
(see Seitz et al. 1998). This is the so called mass-sheet 
degeneracy. Thus in order to get a realistic convergence 
map, A has to be determined by forcing R = at the survey 
scale. The estimate of n is obtain by a non-parametric 
least x 2 method (Bartelmann et al. 1995). The image is 
sampled on a N x N grid, and the potential ona(JV + 2) 2 
grid. Starting from a guess on tp, a guess on the reduced 
shear is obtained and the following function is minimized 
with respect to tp, 



X 2 = ElsW-.9l 



(41) 



The finite difference schemes which are used in order to 
calculate the second derivatives of the potential at pixel 
are, 

Mi, 3) = + 2,j) +?/<(*- 2, j)-2?/>(m)]/4A 2 
Mhj) = [tp{i,j + 2)+tP(i,j-2)-2tP(i,j)]/AA 2 
MhJ) = [^ + l,J + l)-^ + l,.7-l) 

- 1,3 + 1) + 1>(i - 1,3 - 1)]/4A 2 , (42) 

where A is the pixel size. We found that these schemes 
give the best regularization at small scales and avoid the 
usual high frequency oscillations in crude reconstruction 
schemes. At the edges of the field the shift of 2 pixels in 
Eqs. (42) must be only one pixel in the direction perpen- 
dicular to the border since the potential is sampled on a 
(N + 2) 2 grid only. This has the consequence to slightly 
increase the noise at the boundaries but does not produce 
any bias. Once Eq. (41) is minimized, the k map is found, 
and then the condition R = is imposed on it. 



Construction of the initial g map 

A problem that we have to solve in this work is how to 
get the shear pattern of a projected mass distribution? 
Namely from the simulated k we want to get the corre- 
sponding distortion map, put a noise on it, and reconstruct 
k. The construction of a distortion map from a conver- 
gence map is unfortunately an under-constrained problem 
again. For example any transformation of the potential 
tp — > tp + tp with Aip = leaves the convergence un- 
changed, but may change the shear (such a solution is 
for example # 2 — # 2 + 9 x 9 y ). A peculiar solution for the 
potential ip can be obtained by minimizing the \ 2 func- 
tion: 



X 



'"'gUCSS 



(43) 



The resulting reduced shear is then g + g, where g is the 
unphysical solution given by ip, but it is possible to recon- 
struct the convergence using Eq.(41), since k is unchanged 
by the presence of the term g. Unfortunately it is no longer 
the case when the noise is included via Eq.(10), because 
g explicitly comes in the denominator. Fortunately, the 
role of the denominator in Eq.(10) is weak (this is a one 
percent effect on each galaxy if we take g ~ g ~ 0.01 and 
~ 0.1), this is in particular one of the reasons why the 
weak lensing approximation works so well. We thus did 
not try to correct for the presence of a spurious contri- 
bution g in all the calculations, since it gives a negligible 
contribution to the signal (in terms of power spectrum and 
moments). 

Noise generation 

Once the true g maps are obtained the noise is introduced 
in a realistic way, using Eq.10: a sample of background 
galaxies with random intrinsic orientations is sheared, 
from which the k map is reconstructed. The galaxies are 
observed in a grid of superpixels of 2.5 arcmin size (the 
minimum size of the k map simulations), in which the 
number of galaxies JVj per superpixel i is known. In prin- 
ciple this number suffers of the amplification bias (depend- 
ing on the line of sight matter quantity), but this effect 
is neglected here. Ni follows a Gaussian distribution of 
mean N p and variance ^N p . Each pixel i of the image 
gives a local estimate of the reduced shear from the Ni 
galaxies, each of them having an intrinsic ellipticity 
which contributes as a source of noise for the shear signal. 

truncated normalized Gaus- 



3 



sian defined over the range [0, 1], 




p{e) ) oc exp 



where we choose a e = 
typical axis ratio of 0.8). 



(44) 



0.12. (which corresponds to a 
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^From Eq. (10), the observed mean ellipticity e, of a 
number Ni of galaxies in the image plane is given by 



g in the supcrpixel i (Schramm & Kayser 1995, Schnei- 
der & Seitz 1995). A realistic estimate of N p depends on 
the observational context, the telescope used, the opti- 
cal filter, the atmospheric conditions. By choosing N p = 
30 gal/arcmin 2 or 50 gal/arcmin 2 we adopted a conserva- 
tive and reasonable assumption about the telescope time: 
the former is accessible in the I-band with 1.5 hour inte- 
gration at CFHT, while the later is accessible for 4 hours 
in the same conditions. 




(45) 



Where ej is an unbiased estimate of the reduced shear 



This article was processed by the author using Springer- Verlag 
MeX A&A style file L-AA version 3. 



